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We investigate the magnetocaloric properties of certain antiferromagnetic spin systems that have 
already been or very likely can be synthesized as magnetic molecules. It turns out that the special 
geometric frustration which is present in antiferromagnets that consist of corner- sharing triangles 
leads to an enhanced magnetocaloric effect with high cooling rates in the vicinity of the saturation 
field. These findings are compared with the behavior of a simple unfrustrated spin ring as well as 
with the properties of the icosahedron. To our surprise, also for the icosahedron large cooling rates 
can be achieved but due to a different kind of geometric frustration. 
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I. INTRODUCTION 

Antiferromagnetic finite-size spin systems with icosa- 
hedral symmetry constitute very interesting frustrated 
materials with rather unusual magnetic properties. 
Among such properties are jumps to the satura- 
tion magnetization in the cuboctahedron and the 
icosidodecahedroni^ as well as metamagnetic phase tran- 
sitions at zero temperature for instance in the icosa- 
hedron and dodecahedron^^ Some of these proper- 
ties, for instance the large magnetization jump to satu- 
ration, are as well present in the Kagome or other lat- 
tice antiferromagnets. 7 The finite-size systems, neverthe- 
less, have the advantage that due to their smallness many 
properties can be investigated (numerically) exactly with 
possible benefits for our knowledge about infinitely ex- 
tended frustrated spin systems. In this article we in- 
vestigate the magnetocaloric properties of certain spin 
clusters with icosahedral symmetry that turn out to be 
interesting as well. 

Magnetocalorics has a long tradition especially in con- 
nection with cooling by adiabatic demagnetization. The 
first successful attempts to reach the sub-Kelvin region 
date back more than 70 years 4 It is equally well possible 
to extend such an adiabatic process to a full Carnot- 
cycle and thus use the magnetization work for magnetic 
refrigeration, 9 In the past paramagnetic salts have been 
the working medium in both cases which limits the cool- 
ing rate to be not more than 2 K/T. Nowadays gadolin- 
ium compounds such as GdsGasO^ or Gd5Si2Ge2 are 
known to be very efficient refrigerant materials ! 10 : 11 

A unified explanation of a magnetocaloric effect that is 
enhanced compared to a paramagnet is provided by the 
observation that the cooling rate assumes extreme val- 
ues close to configurations with a large excess entropy. 
This can happen at certain phase transitions such as the 
first order [ferromagnetic (I) <-> ferromagnetic(II)] phase 
transition observed in GdsSi2Ge2 compounds?^ at mag- 



netic field driven transitions across a quantum critical 
point) 12 ' 13 or at special values of the magnetic field where 
many ground state Zeeman levels are degeneratei 13 ' 14 ' 15 
In this article we demonstrate that an enhanced 
magnetocaloric effect should be observable in cer- 
tain highly frustrated magnetic molecules of icosa- 
hedral symmetry. We discuss the cuboctahedroni^ 
and the icosidodecahedroni 7 . which are already synthe- 
sized Archimcdian solids as well as the icosahedron 
which in its full symmetry could not yet be achieved 
chemicall y 1 6 1 1 8 1 1 9 Some of the aspects of our inves- 
tigation have been previously discussed in connection 
with the classical^ as well as in the quantum version^ 
of the Kagome lattice antiferromagnet and some one- 
dimensional antifcrromagnets i 13 ' 20 We also like to men- 
tion that magnetocaloric studies have been carried out 
in the field of molecular magnetism recently, but were 
mainly focused on low-field behaviorj^i^^iSii^i 2 ^^^ 8 - 

The article is organized as follows. In Sec. [TT| we dis- 
cuss the basics of magnetocalometry, whereas in Sec. IIIII 
we present the magnetocaloric properties of icosahedral 
bodies and compare them with those of a non-frustrated 
spin ring of the same size. The paper closes with a sum- 
mary. 



II. BASIC THERMODYNAMICS 

A. Heisenberg Model 

The spin systems discussed in this article are modeled 
by an isotropic Heisenberg Hamiltonian augmented with 
a Zeeman term, i.e., 



H = - Juv s{u) ■ s(v) + gn B BS z . (1) 
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s(u) are the individual spin operators at sites u, S = 
y~] u s(u) is the total spin operator, and S z its z- 
component along the homogeneous magnetic field. J uv 
are the matrix elements of the symmetric coupling ma- 
trix. A negative value of J uv corresponds to antiferro- 
magnetic coupling. For the symmetric polytopes dis- 
cussed in the following an antiferromagnetic nearest- 
neighbor exchange of constant size J is assumed. 

B. Magnetocaloric effect 

The magnetocaloric effect consists in cooling or heating 
of a magnetic system in a varying magnetic field. Some 
basic thermodynamics yields the adiabatic (i.e. isen- 
tropic, S = const.) temperature change as function of 
temperature and applied magnetic field, 

(9T\ T fdS\ 

\dBJs ~ ~C(T,B) \dB) T ' 

This rate is also called cooling rate. C(T, B) is the 
temperature- and field-dependent heat capacity of the 
system. For a paramagnet this rate is simply 11 

/ST\ P ™ T 

{ob)s =b (3) 

This situation changes completely for an interacting spin 
system. Depending on the interactions the adiabatic 
cooling rate ^-^ can be smaller or bigger than the para- 
magnetic one and even change sign, i.e. one would ob- 
serve heating during demagnetization and cooling during 
magnetization! 13 i 28 For the purpose of clarity this will be 
shortly illuminated with the help of an antiferromagnet- 
icaliy coupled dimer of two spins with s = 1/2, where 
a singlet constitutes the ground state and a triplet the 
excited state. Following J2]) one notices that in the vicin- 
ity of the magnetic field B c , where the lowest triplet level 
crosses the singlet, the entropy changes drastically at low 
temperatures due to the fact that at the crossing field the 
ground state is degenerate whereas elsewhere it is not. 
This behavior is displayed in Fig. [TJ where below and 
above the crossing field the cooling rate, i.e. the slope of 
the isentropes, assumes large values. 

III. MAGNETOCALORIC EFFECT IN 
MAGNETIC MOLECULES 

Regarding the use as a magnetic refrigerant material 
magnetic molecules possess several advantageous proper- 
ties. They can be synthesized in a great variety of struc- 
tures and they can host various paramagnetic ions. Very 
often they also do not interact magnetically with each 
other in a bulk sample due to large distances between 
the magnetic centers of different molecules that are pro- 
vided by extended ligands, which means that the mag- 
netic properties of a single molecule can be assigned to 
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FIG. 1: (Color online) The curves show isentropes of the an- 
tiferromagnetically coupled dimer of two spins with s — 1/2. 
The straight lines represent isentropes of a paramagnet. Com- 
pared to a paramagnet the cooling rate of the antiferromag- 
netic dimer can be smaller (a), about the same (b), much 
bigger (c), or even negative (d). 



the macroscopic sample. If it would be possible to obtain 
structures with exhibit extraordinarily large ground-state 
degeneracies at certain magnetic fields one could exploit 
these materials for very efficient magnetization cooling. 

Every antiferromagnetically coupled spin system ex- 
hibits ground state level crossings such as the aforemen- 
tioned spin dimer. The cooling rate at such a crossing 
can assume large values, but it turns out that it is pos- 
sible to even increase the rate in a special class of highly 
frustrated magnetic molecules of icosahedral symmetry. 
This will be discussed for the cuboctahedron, which is 
chemically realized as a {Cu^Las} molecule^ and for 
the icosidodecahedron, which is chemically realized as a 
{M072F830} molecule^! and a {M072V30} molecule^ as 
well as for the not yet synthesized icosahedron. The be- 
havior of these frustrated spin systems is compared with 
the behavior of an antiferromagnetic, but not frustrated 
spin ring with TV = 12 sites. All results are obtained by 
means of numerical diagonalization. 

The geometric structures of the discussed bodies can 
be visualized for instance at the following Refs. [I03(J 

The cuboctahedron is one of the smallest antiferromag- 
netic spin systems that can host independent localized 
magnons.— ^ These localized states are intimately con- 
nected with an enhanced degeneracy of energy levels and 
- in extended spin systems such as the Kagome lattice 
antiferromagnet - with the appearance of flat bands. In 
addition, the possibility to arrange several independent 
magnons on the (finite) spin lattice results in a linear 
dependence of the minimal energy E m i n (M) on the total 
magnetic quantum number M. Therefore, at the satura- 
tion field B sat a massive degeneracy of ground state levels 
can be achieved. 

Figure [2] shows that even in a system as small as the 
cuboctahedron the ground state multiplicity at the satu- 
ration field can reach a rather large value, in this case of 
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FIG. 2: (Color online) Low-lying energy levels of the cuboc- 
tahedron with s = 1/2 (dashes) and s = 1 (x-symbols) at the 
saturation field B Bat . The attached numbers give the multi- 
plicity d,M of each M-level. 



nine independent of spin quantum numberpii^ Thus, the 
entropy assumes a non-zero value at zero temperature of 
S = S(T = 0,B = B sat ) = fc B ln(9). All isentropes 
with S < So therefore arrive at the phase space point 
(T = 0,B = B sat ). 
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FIG. 3: (Color online) Isentropes of the cuboctahedron with 
s = 1/2 and s — 1. The entropies here and in the following 
are given in units of Boltzmann's constant fes; from the upper 
left to the lower right the values of S/ks are 0.01, 1.01, 2.01, 
3.01, and 4.01, respectively. B S at = {12s\J\) / (g^s)- 

Figure [3] displays various isentropes of the cubocta- 
hedron both for s = 1/2 (l.h.s.) as well as for s = 1 
(r.h.s.). Both temperature and magnetic field are nor- 
malized to the saturation field here and in the following. 
One clearly sees that the low-entropy curves, S < So, ap- 
proach the .B-axes in a rather universal way independent 
of the magnitude of the spin s. Close to the saturation 
field the slope of the isentropes, i.e. the cooling rate, can 
assume large values. Below the saturation field the isen- 
tropes remain rather flat, i.e. the temperature does not 
increase again when going to B = 0. This is of course an 
important property because otherwise the system would 
heat up again when switching off the field. 

The non-frustrated ring system, Fig. [5J does not pos- 
sess the later property nor does it exhibit large cooling 




FIG. 4: (Color online) Isentropes of the a bipartite, i.e. non- 
frustrated spin ring with N = 12 spins s — 1/2 and s = 1. 
S aat = (8s\J\)/( gf i B ). 



rates at low temperatures close to the saturation field. 
The second property is easily understood because the 
massive degeneracy at the saturation field does not oc- 
cur in an antiferromagnetic ring. Therefore, only isen- 
tropes with S < ks m(2) approach the field axis. The 
first property, however, is related to the overall structure 
of the low-lying levels. In finite bipartite systems the 
levels are arranged in rotational band s 31 ' 32 which means 
that in each sector of total magnetic quantum number M 
the excited states are separated from the ground state in 
the respective sector by a no n- vanishing gap. Thus at 
lower magnetic fields a certain entropy can only be real- 
ized by populating higher-lying levels, i.e. by acquiring 
a higher temperature. 
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FIG. 5: (Color online) Isentropes of the icosahedron with 
s = 1/2 and s = 1. S aat = (14.472s| J\)/(gfi B ). 

The behavior of the icosahedron is intermediate. The 
antiferromagnetic icosahedron is also a geometrically 
frustrated spin system but differs from the cuboctahe- 
dron in some aspects. Each spin has five nearest neigh- 
bors (cuboctahedron - four), and the structure consists 
of edge-sharing triangles instead of corner-sharing trian- 
gles as for the cuboctahedron. Therefore, the icosahe- 
dron does not possess independent one-magnon states 
and consequently a similar degeneracy at the saturation 
field cannot be expected, compare Fig. [5] 

Nevertheless, it turns out that the special frustration 
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FIG. 6: Low- lying energy levels of the icosahedron with s = 
1/2. The attached numbers give the multiplicity cLm of each 
M-level. Note that the energies are rescaled in each sector of 
total spin S. 



of the icosahedron leads to different (quasi) degenera- 
cies in several Hilbert subspaces with total spin S. This 
can be seen in Fig. [H] where the low-lying levels of the 
icosahedron with s = 1/2 are displayed. Note that for 
better visibility the energies are rescaled in each sector 
of total spin 5. The numbers display the multiplicities 
d,M of the M-levels. These are the relevant degeneracies 
in an applied field. The total degeneragy at B = is 
d = c?m • (25+ 1). Thus, we find also for the icosahedron, 
that at certain magnetic field values, not necessarily at 
the saturation field, highly degenerate levels will cross 
and give rise to notable (T = 0)-entropy. This can for 
instance be observed for the isentrope with 5 — 2.01fce, 
which is displayed by the dashed-dotted line in Fig. [5] ft 
approaches the field axis very closely at about 0.9-B sa t- 

The unusually large degeneracy in many sectors of to- 
tal spin is also the reason for the interesting property that 
the isentropes remain at rather small temperatures for 
decreasing magnetic fields. This is in strong contrast to 
spin rings and reflects the fact that populating the degen- 
erate levels produces sufficient entropy without increas- 
ing the temperature. The upturn at small fields close to 
B = is mainly due to the non-degenerate ground state 
in the sector with 5 = which is separated by large gaps 
from excited states. 

The two panels of Fig. [7] summarize the above discus- 
sion for the three systems with N = 12 by comparing 
the isentropes with 5 = 2ks- One sees that at low tem- 
peratures close to the saturation field the cuboctahedron 
indeed achieves the largest cooling rate. In the extreme 
quantum case, i.e. for s = 1/2 it also outperforms the 
two other systems when looking at the achievable tem- 
peratures for B — > 0. Nevertheless, for the icosahedron 
similarly large cooling rates can be achieved due to a 
different kind of geometric frustration. For larger spin 
quantum numbers, i.e. for becoming more classical, the 
differences between the cuboctahedron and the icosahe- 
dron tend to disappear. The unfrustrated ring systems 
always shows poorer cooling. 
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FIG. 7: (Color online) Isentropes 


with S = 2fcs of the cuboc 



tahedron, the icosahedron, and the ring with N = 12 with 
spins s = 1/2 (l.h.s.) and s = 1 (r.h.s.). 
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FIG. 8: (Color online) Isentropes of the icosidodecahedron 
with s = 1/2 and s = 5/2. B sat = (X2s\J\)/{gp B )- 



Finally we like to discuss the behavior of the antiferro- 
magnetic icosidodecahedron which also possess icosahe- 
dral symmetry. This Archimedian solid is closely related 
to the cuboctahedron. It also consists of corner-sharing 
triangles, and each spin has four nearest neighborship 
But compared to the cuboctahedron it has a much big- 
ger ground-state degeneracy of 38 at the saturation field, 
again independent of spin quantum number. Figure [S] 
shows some isentropes of the icosidodecahedron for the 
two experimentally relevant cases of s = 1/2 and s = 5/2 . 
Due to the much larger Hilbert space these isentropes can 
be evaluated exactly only close to the saturation field 
since there only some small subspaces contribute. The 
magnetothermal behavior is very similar to the cuboc- 
tahedron with the noticeable difference that now isen- 
tropes with entropies up to 5o = 5(T = 0, B = B sat ) = 
ks ln(38) ~ 3.63fcs head towards T = at the saturation 
field. 



IV. SUMMARY 

In summary we can say that the investigated frustrated 
antifcrromagnetic bodies show an enhanced cooling rate 
in comparison to non-frustrated (bipartite) spin rings. 
This rate is especially large for those systems that show 
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a large degeneracy of levels, either at the saturation field 
(cuboctahedron) or elsewhere (icosahedron). 

A few words seem to be in order regarding the question 
how realistic the outlined scenario is. In realistic systems 
the perfect degeneracy of levels at the saturation field 
will certainly be lifted, thus keeping the entropy S(T = 
0,-8 = £> sa t) at a small value. Nevertheless, the low- 
energy density of states will remain large in the vicinity of 
the saturation field, since the originally degenerate levels 
move not too far, thus the magnetothermal properties 
will be left qualitatively unchanged, see also the studies 
in Refs. El 



Finally we like to mention that similar effects can be 
observed in interacting electron systems described by the 
Hubbard model. This is related to the appearance of flat 
bands in these systemsi 33 i 34 ' 35 i 36 
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